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INTRODUCTION 


The study of “Evolution” has in recent years been so monopolized 
by experimental geneticists with their interest in induced mutations, 
by cytogeneticists concerned with polyploidy, and by paleontologists 
with their concrete historical evidence, that it may sometimes be 
overlooked that taxonomists have been concerned with phylogenetic 
considerations for at least a century and probably much longer. The 
successive supplanting of Artificial by Natural classifications, and of 
Natural by avowedly Phylogenetic schemes, bears witness to this 
fact. In general, however, the taxonomist treats with such broad 
groups that his attention to any particular population must neces- 
sarily be somewhat perfunctory. To many systematists the classical 
revisionary treatment, with keys, descriptions, and ordered nomen- 
clature, is the ultimate goal of the study. Having fastidiously and 
laboriously completed such a document, they may pass unconcernedly 
on to another, frequently quite unrelated problem. Perhaps the 
taxonomist’s reluctance to suggest or accept premature and too 
facilely constructed genealogies has sometimes led to his dismissal as 
a mere shuffler of specimens and nomenclature, with nothing of value 
to offer to evolutionary theory. 

Although he usually has no evidence from the fossil record and only 
exceptionally has the opportunity to grow and manipulate the plants 
with which he is concerned, the taxonomist may still be able to con- 
tribute something to the evolutionary picture by combining his 
observations on comparative morphology with those gained from 
field study, as it is hoped this paper may indicate. Because he is 


* Introduced by I. M. Johnston and R. C. Rollins. 
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usually concerned with naturally occurring species, he must often 
pioneer, for only recently has the geneticist begun to emerge from his 
laboratory or to shift his attention to other than economically im- 
portant plants. That wild plants of no apparent immediate value 
may reveal additional genetic and evolutionary information, however, 
is quite as likely as that such significant data are discoverable only in 
our cultigens. The taxonomist can render a useful service in pointing 
out and underlining problems suitable for manipulatory and statistical 
attack and he is, at the same time, in a favorable position to test and 
apply hypotheses and generalizations reached on the basis of controlled 
experimentation. It should perhaps be emphasized that when the 
problems he poses are subsequently explored experimentally or mathe- 
matically and the results applied to classification of the organisms 
involved, these activities are as truly taxonomic as his original, more 
exploratory and descriptive studies. Changes in approach and in- 
strumentation do not alter the general objectives of the field. Tax- 
onomy is a sufficiently broad discipline to incorporate contributory 
data from all lines of research, so long as these are properly docu- 
mented and have a bearing upon relationship. 


THE GENERAL PROBLEM IN Cosmanthus 


The present paper is an attempt to distill from a revisionary study 
(Constance, 1949) of Phacelia subgenus Cosmanthus, of the family 
Hydrophyllaceae, some indications as to the evolutionary history of 
three selected entities. An intensive investigation of the group based 
upon comparative morphology, including chromosome number, was 
supplemented by a month’s field work in critical areas, assisted by a 
grant from the Permanent Science Fund of the American Academy of 
Arts and Sciences and the sponsorship of Southern Methodist Uni- 
versity and the Texas Research Foundation. A careful examination 
of dried specimens had made it quite apparent that the interrelation- 
ships of some of the species were too complex to be ascertained by 
herbarium studies alone. The writer hoped that differential chromo- 
some numbers and possibly allopolyploidy might furnish some tenable 
explanations. This hope was not realized because the particular 
entities involved, insofar as cytological material could be obtained, 
all proved to be homoploid and to have nine pairs of chromosomes 
(Cave and Constance, 1947, 1950). 

Nevertheless, the correlation or coincidence of distributions and 
morphological features appears to show a pattern or patterns which 
may prove to be of some general interest. As has been suggested 
above, the revisionary study of this group has posed problems which 
would appear to be suitable for genetic attack, and the study has 
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concluded by showing the eminent desirability of other investigations 
which the writer is, unfortunately, not in a position to undertake. 
It should be firmly borne in mind that the suggestions as to affinity 
made here have no genetic documentation whatsoever and are thus 
necessarily highly circumstantial and conjectural. They may, 
nevertheless, afford some clues to the existence of an evolutionary 
pattern to be sought for in other groups occurring under similar habi- 
tat conditions in the same or completely different areas. 


THREE SPECIFIC CASES IN COSMANTHUS 


Of the three cases selected for presentation here, each involves a 
trio of species or varieties; one member of each trio appears to be 
derived from the other two members of the same trio. These do not 
correspond to the usual situations of two entities of overlapping 
distribution possessing a thin band of sporadic hybrid progeny in the 
zone of overlap, nor to the introgression of characters of one species or 
variety into another one as discussed by Anderson in his recent 
provocative essay (1949) on this subject, although introgression may 
well be involved. Instead, the morphologically intermediate popula- 
tions give every appearance of being as vigorous as their presumed 
parental species, show no cytological irregularities and have distribu- 
tions which indicate that they are successfully established and even ag- 
gressive in areas and habitats occupied by only one or neither of the 
two parental types. If the relationships suggested here are sub- 
stantiated by further studies, they may point the way to modes of 
origin not generally accounted for in the literature of the genetics of 
wild plants. 

No attempt is made to present the relevant data in quantitative 
terms nor to stress their significance by graphic statistical or ideo- 
graphic devices. Statistical techniques, justly or unjustly, often 
carry connotations of objectivity and finality which would be singu- 
larly inappropriate as applied to the purely qualitative data and very 
tentative hypotheses of affinity suggested here. There is no doubt in 
my mind, however, that the employment of hybrid indices of one sort 
or another to the material studied would point in the same direction 
as does the empirical evidence. Instead, it is hoped that a tabular 
presentation of differences and resemblances, accompanied by illus- 
trations of the plants concerned and maps of their distribution, will 
suffice to make the problems clear. 

(1) Phacelia strictifora var. Robbinsii.—Disregarding two of the infra- 
specific populations comprising the notably polytypic Phacelia stricti- 
flora (Engelm. & Gray) A. Gray, we may focus our attention upon 
two others, var. Robbinsii Const. and var. Lundelliana Const., and 
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their relationship to a second species, P. hirsuta Nutt. The characters 
by which var. Robbinsii differs rather strikingly from var. Lundelliana 
are almost exactly those in which it agrees with P. hirsuta (Table 1 


and Figure 1). 


Phacelia strictiflora 
var. Lundelliana 


(n = 9) 


Stems and branches stout. 


Basal leaves persistent, 
succulent, glabrate be- 
neath, entire to shallowly 
toothed. 


Cauline leaves shallowly 
toothed. 


Pedicels strictly erect or 
stiffly ascending, 2-10 
mm. long. 

Calyx lobes, especially 
the lowermost, markedly 


TABLE 1 


Phacelia strictiflora 
var. Robbinsi1 
(n = 9) 
Stems and branches 
slender. 


Basal leaves sub-per- 
sistent, membrana- 
ceous, strigose beneath, 
deeply lobed to pin- 
natifid. 


Cauline leaves deeply 
lobed to pinnatifid. 


Pedicels stiffly to rath- 
er loosely ascending, 
2-10 mm. long. 


Calyx lobes little or not 
at all accrescent. 


Phacelia hirsuta 


(n = 9) 
Stems and branches 
slender. 


Basal leaves early- 
withering, membrana- 
ceous, strigose beneath 
pinnatifid or pinnate. 


Cauline leaves deeply 
lobed to pinnatifid. 


Pedicels loosely ascend 
ing to widely spread- 
ing, 3-15 mm. long. 
Calyx lobes scarcely 
accrescent. 


accrescent. 


Ovules usually 4 to 
each placenta. 


Ovules 8-14 to each pla- Ovules 8-14 to each 
centa. placenta. 
¥ 

In the diagnostically significant features of possession of a basal 
rosette, short pedicels, and a larger number of ovules, it is clear that 
vars. Lundelliana and Roblinsii are members of the same species and 
distinct from P. hirsuta. 

As shown by the map (Figure 2), var. Robbinsii occurs partly near 
but not within the area where the ranges of Phacelia strictiflora and 
P. hirsuta overlap. The habitat of var. Lundelliana appears always 
to be in sandy soil, often in or near oak woods; P. hirsuta has been 
collected chiefly in oak woods, but usually in loamy alluvial, calcareous 
or cherty soil. The habitat of var. Robbinsii is usually noted as the 
loam of open meadows, forest openings, or rocky hillsides; the type 
locality contains weathered material from nearby granite knobs. The 
habitat of var. Robbinsi thus approaches in some respects that of P. 
hirsuta. The areas comprising these habitats have been so greatly 
altered by clearing, grazing, and agriculture that it is not clear what 
the original ecological preferences of these entities may have been. 
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Inasmuch as var. Robbinsw appears to be fertile and self-perpetuating, 
and has extended its range beyond that of both the forms which 
circumstantial evidence suggests may have been its parents, it may be 
supposed that there was some available intermediate habitat (original 
or man-made) in which the hybrid derivatives were more successful 
than were the original forms themselves. 
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Figure 2. Distribution of Phacelia strictiflora vars. Lundelliana and 
Robbinsii and P. hirsuta. 


(2) Phacelia patuliflora var. patuliflora——The second case involves 
Phacelia patuliflora (Engelm. & Gray) A. Gray and P. laxa Small. 
The former species has two phases, var. teucriifolia (Johnst.) Const., 
chiefly of the Edwards Plateau of Texas but extending southward into 
Coahuila, and var. patuliflora of central and southeastern Texas and 
adjacent Mexico. Morphologically, var. patuliflora is extremely 
variable, and appears to consist of numerous interlocking geographical 
races and local forms, whereas var. teucriifolia and P. laxa are both 
relatively uniform. Moreover, if the broad gamut of variations and 
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combinations exhibited by var. patuliflora could be arranged in a 
linear series, P. laxa would stand at one end of the series and var. 
teucritfolia at the other (Table 2 and Figure 3). 


Phacelia patuliflora 
var. teucritfolia 
(n = 9) 
Branches rather stiffly 
ascending. 


Pubescence of stems 
finely and densely strigu- 
lose. 


Pedicels ascending. 


Calyx lobes mostly lan- 
ceolate, acute, ascending 


TABLE 2 


Phacelia patuliflora 
var. patuliflora 


(n = 9) 


Branches usually more 
or less decumbent. 


Pubescence of stems 
finely strigulose to 
densely spreading- 
hirsute. 


Pedicels ascending to 
widely spreading or 
reflexed. 


Calyx lobes narrowly 
oblong or lanceolate to 


Phacelia laxa 


(n = 9) 
Branches diffuse, pros- 
trate to ascending. 
Pubescence of stems 
sparsely spreading- 
hirsute. 


Pedicels widely spread- 
ing or reflexed. 


Calyx lobes oblanceo- 
late to obovate, obtuse, 





rotate-spreading to re- 
flexed in fruit. 


in fruit. obovate, obtuse or 
acute, spreading to as- 


cending in fruit. 


Corolla small, broadly 
campanulate. 


Corolla variable in 
size, rotate-campanu- 
late to broadly cam- 
panulate. 


Corolla large, rotate- 
campanulate. 


Some of the populations of var. patuliflora, particularly in the Rio 
Grande Plains, so nearly simulate var. teucriifolia (perhaps as a 
consequence of introgression) that it seems unwise to attempt to 
separate them specifically. On the other hand, it is usually possible 
to achieve a satisfactory separation from the wholly sympatric P. laza, 
at least in the field. 

In this instance the ranges of the two presumed parents do not 
intersect, at least at present, but they are broadly connected by that 
of var. patuliflora, which is postulated to comprise the sum of their 
existing hybrid derivatives. The usual habitat of var. teucriifolia 
appears to be alluvial soil along streams and in draws in West Texas; 
P. laxa is characteristically found in moist alluvial thickets along the 
lower courses of rivers in the Nueces Bay region (Figure 4). The 
var. patuliflora is most common in sandy soil, especially on river 
terraces and coastal dunes, but it occurs over a wide variety of habi- 
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tats, including roadsides and railway embankments, as well as alluvial 
thickets and open, sandy woods. It seems especially probable in this 
case that man’s destruction of old habitats and creation of new en- 
vironments has permitted establishment and spread of a new entity. 
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Figure 4. Distribution of Phacelia patuliflora vars. teucriifolia and 
patulifiora and P. lara. 


(3) Phacelia gilioides.—The third selection involves Phacelia gilroides 
Brand and its affinities with P. Purshii Michx. and P. hirsuta Nutt. 
The species of the genus occurring in Missouri have been confused 
because attempts were made to treat these three entities as only two. 
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Very different ideas of the distributional patterns of P. Purshii and 
P. hirsuta were thus obtained, depending upon which morphological 
criteria were employed for separation. The recognition of a third 
species in the area greatly assists in clarifying the picture. This 
entity, P. gilioides, is suggested to have arisen as a consequence of 
hybridization between the other two, but it appears now to be a 
fully fertile, self-perpetuating, and morphologically stable species, 
which has colonized some areas uninhabited by either of its hypotheti- 
cal parents. The morphological evidence for this suggestion as to 
relationship is summarized in Table 3 and Figure 5 and the geographi- 


cal in Figure 6. 


Phacelia hirsuta 
(n = 9) 
Pubescence of stems and 
inflorescence spreading- 
hirsute. 


Lobes of the cauline 
leaves mostly obtuse. 


Corolla lobes entire, pi- 
lose on the back. 


Ovules usually 4 to each 


TABLE 3 
Phacelia gilioides 
(n = 9) 
Pubescence of stems 
and inflorescence stri- 
gose, the latter canes- 
cent. 
Lobes of the cauline 
leaves usually acute. 
Corolla lobes denticu- 
late to fimbriate, pilose 
on the back. 
Ovules usually 4 to 


Phacelia Purshii 
(n = 9) 

Pubescence of stems 
and inflorescence stri- 
gose or strigulose, often 
sparse. 
Lobes of the cauline 
leaves acute. 
Corolla lobes deeply 
fimbriate, usually gla- 
brous. 
Ovules only 2 to each 
placenta. 





placenta. each placenta. 


The writer has no first-hand knowledge of ecological conditions 
where these three species occur in Missouri, and herbarium labels are 
not especially revealing. Phacelia Purshii is apparently confined 
largely to alluvial bottomlands; P. hirsuta appears to be abundant in 
woodlands, especially in cherty soil. The ecological niche most 
commonly assigned to specimens of P. gilioides is limestone glades or 
“barrens,” but whether these sites differ appreciably from those 
supporting the partially sympatric P. hirsuta is not clear, and the two 
have been mixed, although rarely, in the same collection. Again, the 
ranges of the supposed parents are not known to overlap, but that of 
the assumed hybrid derivative connects them. Field study might 
reasonably be expected to provide some clue as to the environmental 
reasons for the persistence and successful spread of P. gilioides. 


DISCUSSION AND CONCLUSION 


The linear arrangement of species adopted by the writer in his 
recently published revision of the subgenus Cosmanthus had no special 
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significance. However, it was an attempt to approximate an assumed 
phyletic progression, giving particular weight to reduction in number 
of ovules—which appears to be a general trend in the genus—and to 
suspected migrational history. If this arrangement has any merit, it 
tends to show that the species may be grouped into about four levels 
of advancement. Of those with which we are presently concerned, P. 
strictiflora and P. patuliflora would occupy the second level, P. hirsuta 
and P. laxa the third, and P. Purshi the fourth. This observation is 
repeated here only because it indicates that the supposed cases of 
hybrid derivation discussed above would all necessitate the crossing 
of comparatively distantly related parents. It would, perhaps, be 
less surprising if, for example, P. patuliflora crossed with P. stricti- 
flora, which is of about the same level of advancement and overlaps 
it broadly in distribution, but this does not appear to occur. The 
hybridization which is significant in this group, in terms of the 
success of the resultant offspring, appears rather to be between species 
which are well defined morphologically and not too closely related. 

These three examples will, it is hoped, indicate that the taxonomist, 
working only with his ordinary tools, may be able to throw some light 
on the phylogeny of the groups with which he is concerned. In none 
of the cases cited has the proposed explanation been established, and 
supplementary genetic or other investigation is earnestly invited. 
Only by evidence from every available source, guided by the fitful 
illumination of successive and even contradictory hypotheses, can we 
hope to approach a truly phylogenetic classification. 
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SUMMARY 

This is an attempt toward a comprehensive and rigorous treatment 
of acoustic radiation pressure in the case of non-viscous motions. 
Unlike the optical radiation pressure, for which exact formulas can be 
derived, the acoustic radiation pressure is usually computed by 
approximation methods, the accuracy of which depends on the 
“amplitude” of the motion. The theory presented in this paper 
begins with a rigorous definition of the amplitude of a whole motion, 
in contrast to the amplitude of vibration at a particular point. From 
there on the first- and second-order terms in the series expansion of 
the acoustic radiation pressure, and other quantities, with respect to 
the amplitude are rigorously derived. The development of the theory 
brings out a remarkable mass-energy relation: defining the equilib- 
rium of a volume | of the fluid by the average position of its bound- 
aries (with respect to time), the difference (, — E;) of the potential 
and kinetic energies inside this volume V is equal to 4 mc*, where c is 
the wave velocity and m the time-average increase of mass caused by 
the motion inside the volume V’. Finally, a formula is obtained for 
the radiation pressure on the walls of a container which participates 
in the oscillations of the fluid. Some of the formulas involve a higher 
approximation than that obtainable from the wave equation V’¢ = 
¥/c? (g = velocity potential), and call for a more rigorous solution of 
hydrodynamic problems concerning radiation pressure. The classical 
expressions for the force on a completely immersed body are not 
affected, however, as they only depend on the differences in pressure 
between various points of fluid, and not on absolute values. Ex- 
periments to measure these absolute values are suggested. 


Definitions 


1. The acoustic radiation pressure Ap; at a particular point P, 


+t This work was completed under contract N6onr-262 between the Office 
of Naval Research and Wesleyan University. The author is now at The 
Catholic University of America, Washington, D. C. 

* Introduced by W. G. Cady. 
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fixed in space will be defined as the average change in pressure (p;,. — 
pie) from the reference pressure p;,o at this point: 








ADs = Pit — Dio (1) 


This definition is useful only when the average is practically inde- 
pendent of the length of time over which it is computed, as, for in- 
stance, in the case of continuous oscillatory motions (and also steady 
flows). Its physical meaning can be visualized by imagining that a 
small area dS of the walls (assumed, for simplicity, to be rigid and 
stationary) containing the fluid is replaced by a piston of large inertia. 
Then a constant external force dF has to be applied to this piston to 
keep it from moving away from the wall. According to the law of 
conservation of momentum, this force is given by the relation: 


dF = Ap; X dS (2) 


which shows that the radiation pressure is equal to the force acting on 
the unit area of the walls. 

Experimentally the acoustic radiation pressure is always measured 
by its action on a wall or solid surface, and therefore is often computed 
by using expressions that are correct on the walls, but not in the bulk 
of the fluid. Furthermore, when the walls are those of a completely 
immersed body, the total force acting on this body according to Eq. 
(2) will not be changed by adding a constant to any expression for the 
radiation pressure, which again allows more freedom in the choice of a 
particular formula. 

2. The average increase in density Ap; at a point P; fixed in space 
is defined in the same way as the average increase in pressure, accord- 
ing to the relation: 





Api = pis — Pio (3) 


This function is introduced here chiefly as a means of expressing mathe- 
matical relations in a more convenient way, its experimental signifi- 
cance being less obvious than that of the acoustic radiation pressure. 
In the case of a transparent fluid enclosed by walls in which two plane 
parallel glass windows are provided, however, photometric measure- 
ments using optical interference or absorption might make it possible 
to measure directly the average of the above function Ap; over the 
volume of a cylinder having the two: windows for bases. This would 
make it possible to check experimentally Eq. (33) below, which in- 
volves the average increase of mass caused by the motion in such a 
volume. 
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3. In the present treatment the displacement £,, (t) of any particie 
m at time ¢t is regarded as a function of an amplitude A, but not neces- 
sarily as proportional to A. At any given value of A, £,, (¢) is com- 
pletely determined when the nature and configuration of the system 
are known. This particular value of A defines one member of an 
infinite family of motions. 

The amplitude A is a parameter that may be defined in any one of 
a great variety of ways. It may be the amplitude of motion of a 
piston, or the amplitude of the excess pressure at a given point, or the 
amplitude of the emf that drives an electromechanical transducer; on 
the other hand, we prefer to regard it as a generalized amplitude, and 
not necessarily as any physical amplitude, as long as it can be used in 
a power series of the type shown in Eq. (4) below. That which 
concerns us here is not the mechanism by which the acoustic waves 
are produced, but their variation with amplitude. 

The configuration of the system is assumed to be completely 
specified by a set of coefficients &, (m,t), where h = 1, 2, 3, . . . 
Each coefficient is a function of time for any particle m, and the aggre- 
gate, for all values of m and h, describes the entire system, including 
the boundary conditions. The &, are therefore characteristic of the 
family, and are the same for all members, that is, for all values of A. 

We now express the displacement £,, (¢) of any particle m as a 
convergent power series in A: 


En(t) = A&, (m, t) + A*k (m,t) + . . + AME (mt) +... (4) 


in which ¢ is the time and m the symbol of the particle, usually specified 
in terms of its coordinates 2m, Ym, 2m in the equilibrium state. 

The method indicated above systematizes a procedure implicitly 
contained in all treatments of wave problems. It is common practice 
to include in the expression for the displacement a factor representing 
the physical amplitude of the motion at some point, such as the 
source of the waves. In our definition, each value of this physical 
amplitude at the source corresponds to and characterizes a different 
state of motion. The physicist accustomed to the first-order approx- 
imation might contend that changing the physical amplitude at the 
source only results in multiplying the displacements at all other 
points by a common factor, and that the difference is entirely trivial. 
To this the answer is that we are dealing here with a more accurate 
picture of the motion. Increasing sufficiently the physical amplitude 
at the source will cause a distortion represented by the terms in A?, 
A’, etc. in Eq. (4), and therefore the new state of motion will be very 
different from the old. Moreover, our definition does not require the 
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generalized amplitude A to coincide with any physical amplitude, as 
long as it can be used in a series expansion of the type illustrated in 
Eq. (4). 

4. When the fundamental laws of mechanics are applied to the 
fluid, the resulting equations can be written as power series of the 
generalized amplitude A equated to zero. An exact solution would 
make each term of these series equal to zero. The theory of acoustic 
radiation pressure, however, only deals with the terms of the first and 
second degree with respect to the amplitude. Even so, each result 
has to be carefully scrutinized to make sure that no second-order term 
is omitted in its derivation. 

5. This use of series expansions in terms of a generalized amplitude 
A can afford a most straightforward and powerful means of solving 
acoustic problems. As an illustration, let us consider the case of 
sound waves in tubes and horns, the solution of which is well known 
because of its importance in the construction of loudspeakers.** We 
start from Lagrange’s law of motion: 


" 7) re) 
pk = — p’ (dot = —, prime = —), § 
é J ot 4 on? (5) 


in which 2 is the distance from the source of the waves, and from the 
compressibility law: 


P — Po = — poco’ [t’ + DEF +... «| (6) 


in which D is what we may call a “distortion coefficient,” since it 
would be zero for a fluid obeying Hooke’s law. On replacing & by 
its series expansion in Eq. (4), and p by its expansion in Eq. (6), 
Lagrange’s equation splits into the new relations: 


é, = i Se (7) 
bo = co? [&”” — D (&2)'] (8) 


The solution of the first equation is the undistorted, fundamental 
wave representing the first order approximation. Assuming this 
fundamental wave to be simple harmonic we replace & by e## ‘¢~ 2/0 
in the second relation, thus obtaining, for the amplitude of the double- 
frequency harmonic in the complete motion, a differential equation of 
the type 


fy)" + f(y) = (9) 


with y = 2waz/c. This equation is integrated by using the particular 
solution f(y) = ye~” , which yields at once the correct expression for 
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the double-frequency harmonic generated by the distortion of the 
fundamental wave: 


» 


i“ 


s Ah = — "Dx etia(t-z/e) (10) 


In this equation A is the displacement amplitude of the fundamental 
wave, and has been used throughout as generalized amplitude for the 


complete motion. 
6. The quantities that have previously been defined, 2. e., the 


acoustic radiation pressure Ap; and the average increase in density 
Ap;, are governed by two fundamental relations: the first deals with 
the static properties of the fluid and can be derived from its law of 
compressibility alone; the second is kinetic in nature and is obtained 
when the fundamental laws of dynamics are applied to the motion. 
These relations will be treated in the following paragraphs. 
Static Relation 

7. The law of compressibility of the fluid can be expanded as a 
power series of Ap = (p — po), in which pp is the pressure in the 
reference state: 


p=pth(Ap)+h(Ap’?+... (11) 
If we introduce the wave velocity, 
dp 
dp 


(12) 


and the distortion coefficient D already occurring in Eq. (6), this 
expression of the compressibility law becomes: 


l D— 


l 
Ap — (Ap)?+... (13) 








Pp = pot “ 


poc* 


8. The distortion coefficient D is zero for a fluid obeying Hooke’s 
law, 2. e., when the pressure is a linear function of the volume. In the 
case of sound waves in gases, if the compression is adiabatic, the value 
of the distortion coefficient D is (y + 1)/2, in which y is the ratio of 
the specific heats. This gives D = 1.20 in the case of air. In 
liquids the distortion coefficient may be much larger; D = 5 in the 
case of water. As will be shown later (Par. 15), however, the larger 
displacement amplitudes experimentally obtainable in the case of 
gases more than compensate for the smaller distortion coefficient in 
producing observable effects. 
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9. Taking averages in Eq. (13) gives the new relation: 


_ _. «-«‘D-1 
Ap = c’*Ap + 





(Ap)? (14) 
Po 


Observing that poc? is the compressional modulus / of the fluid, and 
introducing the average potential energy per unit volume, 


(Ap)? 
= (15) 

we can write Eq. (14) as: 
Ap = “Ap + 2(D — 1) @, (16) 


which is the “static relation’’ between the acoustic radiation pressure 





Ap and the average change in density Ap. 


Kinetic Relation 


10. The fundamental laws of dynamics, applied to a fluid, may be 
expressed in the form of Euler’s equation. For the present purpose 
this equation is preferably written in vectorial form: 


dp Ov v | 
H-d—-—= (= + grad — — ¢ X curl ‘) - (ll (17) 
: p ot 2 

The theory of the acoustic radiation pressure supposes the external 
field H negligible, so that the first term on the left-hand side can be 
omitted. The last term on the right-hand side, v X rot v, is of the 
third degree with respect to the amplitude 4, since the curl of 
the velocity is always a second-order quantity. This last term, 
therefore, is not to be considered. On the left-hand side 1/p can be 
expanded into a power series as in Eq. (13), which gives: 


] ] ] ; 
—=—|1- " (p — Po) + —— | (18) 
p Po 





Po 
and, since the field H/ has been assumed to be negligible, and therefore 
po constant over the fluid, 


di 1 (p — po)’ | 
hr = hb scicu (19) 
p Po poc™ 





When the average with respect to the time is taken, Eq. (19) becomes: 


‘ini ] — 
(dp/p) = —d(Ap—@&+...) (20) 
Po 
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in which é, is the potential energy defined by Eq. (15). The term 
dv/dt on the right-hand side of Eq. (17) has a zero average since the 
velocity does not increase indefinitely. Consequently, if we also 
introduce the average kinetic energy per unit volume 


x = Vo por, (21) 
it follows that Eq. (17) can be written 
— d (Ap — @,) = da, (22) 


and an integration gives the “dynamic relation’”’ 
Ap =C+é6,—-— & 23) 


in which C has the same value at all points of the space occupied by the 
fluid. This constant C disappears when the difference in pressure 
between two points of the fluid is computed, and as a result Eq. (23) 
has been widely used to compute the forces acting on solid bodies 
completely immersed in a vibrating fluid.44; It can also be applied 
to the case of a steady-flow motion, and leads then directly to Ber- 
noulli’s theorem. 


Absolute Value of the Acoustic Radiation Pressure 


11. By using the static relation (16) the value of the constant C in 
the dynamic relation (23) can be computed from the boundary con- 


ditions of the fluid. To that effect the radiation pressure Ap is first 
eliminated between the two relations. This gives: 


C = Ap — (3 — 2D) & + & (24) 
Let us now introduce an arbitrary volume V limited by a surface S 
fixed in space and completely immersed in the fluid. When Eq. (24) 
is integrated over the volume, the integration of Ap will give the 
average mass increase AM caused by the motion in this volume V. 
Denoting the total potential and kinetic energies by FE, and FE; we have: 


br cane vad 
Ca— eAM — (3 — 2D) E,+ Bs | (25) 


12. The average mass increase AM in Eq. (25) can be computed 
from the conditions at the boundary S of the volume VV’. The parti- 
cles which lie on S in the reference state of the fluid are on a surface S’ 
at the instant ¢. Let h be the distance from a point P of S to the 
point P’ of S’ situated on the same normal to S, the positive sense on 
this normal being toward the outside of the volume JV" (see Fig. 1). 
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Let p be the density at the same point P on the surface S. Then the 
average mass increase AWM is given by the approximation formula 


AM = — f phdS (26) 
S 


which applies exactly to all terms of the first and second degree with 
respect to the amplitude A, but not to terms of higher degree. 





13. The calculations leading to Eq. (26) illustrate the statement 
made in Par. 4 about the necessity of discussing carefully the degree 
of approximation achieved in each case. Referring to Fig. 1 again, 
let zh be the distance from an element @ of the cylinder of base dS 
and height h to the surface S, x being a dimensionless number less 
than unity. Then the mass contained in the cylinder is 


1 
dM = aS f pghi (27) 
0 
] 
= hdS f podx (28) 
0 
We can expand pg as: 
Op 
pg=p+—het+... (29) 
On | 


in which p is the density at the point P and dp/dn the derivative of this 
function p along the normal to S. Now p can be expanded in terms of 
the amplitude A, according to the formula: 
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p=ptp(t,y,z0At... (30) 


in which pp is independent of x, y, z, and ¢ since the external field is 
negligible as stated in Par. 10. Therefore dp/dn will start with a 
term of the first degree in A and we have 


Op) 
On 
By use of this relation, Eq. (28) yields, after integrating, 
Opi ' 
dM = (1 + Ph a A+... ) dS (32) 
n 


From Eq. (4) it follows that the series expansion of h starts with a 
term of the first degree. Therefore the second term in the parenthesis 
is of the third degree in A and can be neglected. The mass of the 
elementary cylinder can therefore be written as merely equal to hp, 
which leads directly to Eq. (26), giving the average mass increase 


AM in the volume V. 

NOTE: Owing to a compensation which develops when the inte- 
gration is carried out over the whole surface S it is permissible to use 
the normal component &, of the displacement of the fluid particles, or 
the normal component (&), of the displacement of each point of the 
boundary, instead of the normal distance A in all relations starting 
from Eq. (26). 

Practical Applications 

14. Summarizing Eqs. (23), (25), and (26) the acoustic radiation 

pressure can be expressed by the formula: 


— 1 a - ™ 


AM = — f phdS (34) 
S 


with 





In this equation the various energies E,, Ex, ep, ex, when terms of 
degree higher than the second are disregarded, depend only on the 
first order terms in the expansion of the displacement, according to 
Eq. (4). Therefore all these energies can be computed from the well- 
known wave equation which governs the first-order terms in the 
series expansion of the displacement. This wave equation is best 
written by using the velocity potential y, and then takes the form: 
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Cop fp Hey 1 eo 


(2 or 


(35) 


eee le 


022 Oy?—s 2? 


The situation is different, however, in the case of the mass increase 


AM appearing in the first term of Eq. (33). The displacement h of 
the boundary is of the first order in A but the density contains a 
constant term po, and we can write Eq. (34) in the form: 


——— 


aM =— f wtmdt. . .)(mA+ beA? +...) dS (36) 


From this relation it is apparent that a first-order approximation will 
be satisfactory as far as the density p is concerned; but since the 
second-order term in A combines with po, this term must be known 
exactly if the desired second-order approximation is to be obtained 


for AM and, consequently, for the radiation pressure Ap. The wave 
equation (35), therefore, is not in general sufficiently accurate to be 
used in computing the absolute value of the acoustic radiation pressure. 

15. There are cases, however, in which it is not necessary to solve 
any wave equation in order to know the displacement h of the bound- 
aries. One of these is when the boundaries are made of some solid 
material subjected artificially to a motion practically independent of 
the reaction of the fluid: a sounding pipe closed by a piston to which a 
heavy balance-wheel imposes a well defined oscillatory motion offers 
such a situation. In that case h is exactly known, and solving the 
first-order wave equation (35) will be sufficient to obtain the first- 
order terms needed in Eq. (36) to express the density p. The motion 
of a fluid enclosed within rigid, stationary walls also belongs in the 
same class: in that case h is identically zero and therefore AM is zero. 
Only the energy terms are left in Eq. (33) for the radiation pressure, 
which becomes: 


onal ] . 
Ap = Vy |e ne 3) E, + E, + €p ~~ Gp (37) 


NOTE: In the particular case of plane simple-harmonic standing 
waves filling entirely a rigid tube closed at both ends, the quotients 
E,/V and E;/V are both equal to 1% (é;) max, since the distribution of 
energy is sinusoidal. At a node of pressure, é, is then zero while é, 
is a maximum. Equation (37) then yields for the radiation pressure 
at the nodes: 


D — 2 (Apmax)’ 





Ap —_ (D = 2) (x) max —_ (D —_ 2) (i Dscen _ 4 ke (38) 
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If a small orifice is made in the walls of the tube at this node of pres- 
sure, a certain volume V of fluid will go in or out according to the sign 


of Ap, and we have the relation: 


AV Ap 2-—D /Apmax\? 2-—D 

a ( ) —_ (Ouax)” (39) 
V k 4 k 4 

in which @ is the dilatation. In the case of a gas the maximum dila- 

tation 6.x can be made fairly large and it might be feasible to test 

this relation experimentally. 

16. Another case in which the first-order approximation given by 
the wave equation (35) is sufficient to calculate the acoustic radiation 
pressure is when an acoustic experiment is made in a very large space 
with absorbing walls so that the only important part of the phenome- 
non is in the relatively small portion of the space in which the instru- 
ments are located. ‘Then Eq. (33) reduces to 


Ap = lp — Ck (40) 


This can be seen by referring to Eqs. (23) and (25): The constant C in 
Eq. (23) is given by Eq. (25) and can be written 


C= c?(Ap) av. — (3 — 2D) (@p) av. + (€z)av. (41) 


the subscript “‘av.”’ meaning average over the volume V. In a very 
large space all these averages are zero as long as the waves are not 
reverberated into some troublesome standing pattern. In this 
respect it should be observed that absorbing walls can always be re- 
placed by solid surfaces subjected, artifically if need be, to a suitable 
motion, so that the theory developed here automatically applies. 








The Mass-Energy Relation 


17. Expression (33) for the acoustic radiation pressure can be 
written in another form in which all terms can be calculated from the 
first-order approximation represented by the wave equation (34). 
This is achieved by selecting an equilibrium state of suitable proper- 
ties. In Eq. (4) for the displacement there is no assumption concern- 
ing the equilibrium state. In all mechanical problems the equilibrium 
state is primarily a reference state and is selected merely for con- 
venience. In the present case it need not be the same for all the 
motions of the family (that is, for all values of the amplitude A). 
Not only the equilibrium density pp may vary, but also the shape of 
the container and therefore the relative positions of the particles. 
This being so, let us define a “fundamental equilibrium state” by the 
condition that the average value of the displacement h of the bound- 
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aries during the motion, as defined in Par. 14, be zero. In other 
words, let the reference state be defined by the average position of the 
boundaries. Then it can be shown (see Appendix) that the average 
mass increase m, with respect to the fundamental state, is given by 
the mass-energy relation: 


lg me? = E, — E, (42) 


Using this expression in Eq. (33) we obtain the following expression 
for the acoustic radiation pressure, with respect to this fundamental 
state: 


p=1/V (2D—-1) E£, — EE.) +4 - & (43) 


This involves only the potential and kinetic energies and can therefore 
be calculated from the wave equation (35) for the first-order terms of 
the displacement. This is not in disagreement with the conclusions of 
Par. 14, but merely a different way of expressing these conclusions. 
It is now necessary to find the average position of the boundaries to a 
second-order approximation, and a discussion of this problem leads 
to the same results as stated in Pars. 14 and 15. 

18. The derivation of the mass-energy formula (42) is given in the 
appendix. This formula can also be checked in cases simple enough 
for a complete theoretical treatment: the simplest case is that of a 
rigid container communicating with a cylinder in which a piston is 
subjected to a motion that is infinitely slow, but otherwise arbitrary; 
in that case there is no kinetic energy and the potential energy per 
unit volume is the same at all points. We have also treated the most 
general case of plane waves in a sound pipe, the case of progressive 
spherical waves, and finally, several cases of standing spherical waves. 
In all these cases the mass increase m was computed from the displace- 
ment and the dilatation at the boundaries of the fluid, using Eq. (26). 
These computations are interesting but too lengthy to be reproduced 
here. 

Order of Magnitude of the Acoustic Radiation Pressure 


19. The variation of the acoustic radiation pressure from one point 
to another is always a second-order quantity, according to Eq. (23). 
The absolute values, however, according to Eqs. (33) and (36), are in 
general of the first-order with respect to the amplitude A. They 
become of the second-order if the equilibrium state is defined by the 
average position of the boundaries, to the first-order approximation 
given by the wave equation (35), since the coefficient h; is then zero 
in Eq. (36). If this same definition of the equilibrium state is carried 
out to a higher degree of approximation, then the acoustic radiation 
pressure is not only a quantity of the second-order but the quantity 
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defined by Eq. (43). It can be shown that when the equilibrium state 
is changed so that the equilibrium pressure changes from pp to po + 
6po, the formulas derived in this paper all lead to a radiation pres- 


sure Ap — dp, as should be expected. 


Average Pressure on a Moving Boundary 


20. The acoustic radiation pressure Ap has been defined as the 
average change in pressure at a point fixed in space. If & is the dis- 


placement of an element of the boundary, and Ap, the change in 
pressure on this element, with respect to the reference pressure, we 
have: 
a  achis Ov | 
Ap, = Ap+ grad p - & = Ap— p = Ep (44) 





in which Ap is the change in pressure at the point of reference. Inte- 
grating by parts gives 


t 


ae ] 
lite ta p> & + pv - Bw (45) 








The second term on the right-hand side is zero since both v and & 
are limited quantities, and therefore we are left with: 


Ap = Ap+ pvt (46) 


When the motions of the walls and of the fluid are identical, the cor- 
rective term is twice the kinetic energy é, per unit volume. If either 
v or t% is normal to the surface of the boundary, the scalar product 
» - % 1s the square of the normal velocity. If the boundary is sta- 
tionary, there is no correction to be made, no matter what the tan- 
gential velocity of the fluid may be. 

21. As an illustration, let us consider the case of a single progressive 
plane wave traveling in a pipe. This case can be idealized by imag- 
ining two oscillating pistons, one at each end of the pipe, with motions 
identical to those of the progressive wave. Since the potential and 
kinetic energies are equal in this type of wave, and have the same 
value at all points, Eq. (43) reduces to: 


Ap =2(D—1)& (47) 


which is the value of the average pressure on the stationary walls. 
From this the average pressure Ap, on the pistons is calculated by 
using Eq. (46), which gives: 


Ap, = 2Dé, (48) 
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These results can easily be checked by using the complete solution 
already described in Paragraph 5. 
APPENDIX 


The mass-energy relation (42) can be derived by starting from the 
expression for the kinetic energy E, in terms of the velocity potential 


l d¢\? d¢\? d¢\? 
ot [FEO « 
2 Ox Oy Oz 
v 
in which V is the volume enclosed by the surface S, and dV any 
element of this volume. This expression can be written 








hee tier we oy?e dV (50) 
4b = ——" Do (7 — — oe ae 
os (+ aI ) “‘. ev’? \ 
v v 
in which dl stands for any linear element. Converting the first 


integral into a surface integral and transforming the second by using 
the wave equation (35), we obtain 


1 dy 1 po — By 
eee ine 
ete ie 22) © a 


dV=h—f, (51) 





An integration by parts gives 
1 iH dg |' 1 po d9\? | 
cone Pn a (=) dV — (52) 


he wien 
7 2 @ t 
V 





In the case of an oscillatory motion the velocity potential is a limited 
function of the time ft, since the equilibrium state has been defined by 
the average position of the particles. Therefore the first integral is 
zero. The second represents the potential energy of the fluid, so 

that Eq. (51) becomes, after suitable transpositions of terms, 
E,-—E = ds 53. 
~*~ —— fh e— ds (53) 

° 2 on 
JS 

In this relation dg/dn represents the normal component of the velocity 
at any fixed, geometrical point of the surface S. But it can be shown 
in general that such a function of a geometrical point differs from the 
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corresponding particle function only by quantities of orders higher 
than its own. Therefore, continuing to neglect all terms of order 
higher than the second, we have 
_— = 
_ ee eee g 
Aiea 2°" at 
Ss 





dS (54) 


in which £, denotes the distance traveled by the particle, normally to 
the surface S, from its reference position. Integrating by parts gives 
a zero term again as in relation (52) and we are left with 


1 0¢ = 
E, — E; 3 Po i. dS (55) 
s 
Now 0¢/dt can be calculated from the change in density (p — po) 
since according to the wave equation (35) the first-order approxima- 
tion for the dilatation of the fluid is (p9 — p)/po = (d¢/dt)/c?. Eq. 
(55) can therefore be written 





E, — E, - 4b ¢ (p. — p) En dS (56) 
S 


So far the displacement &, and the density increase (p — pg) in this 
equation are not exactly the actual ones, but only the first-order 
approximation functions satisfying the wave equation Y7¢ = (1/c?)d?¢/ 
df. But the difference between the actual and approximate values 
is of the second-order with respect to the amplitude, and using the 
actual values in Eq. (56) will therefore result in an error of the third 
order. Finally, if the equilibrium state is defined by the average 
position of the walls to a second-order approximation the quantity 


posn Will be zero to this approximation, and the result is: 





1 _— 1 
een a - 
E,-— Ek = -“ c pt, dS = : me (57) 
s 
in which m is the average mass increase in the volume !’ occupied by 
the fluid in the equilibrium state. 
NOTE: We could also write in general, using, Eq. (56), 


l 1 i 2 1 ‘an 
5 CAM a he pt, dS = E, — E, — > Pee En dS (58) 


S S 
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This would make it immediately evident that &, has to be known to a 
second-order approximation if AM is to be computed to this approxi- 
mation. If the equilibrium state is defined by the average position 
of the walls, the last integral drops out and the result is again Eq. (57). 
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THERMODYNAMIC RELATIONS IN 
n-VARIABLE SYSTEMS IN JACOBIAN FORM: 


PART I, GENERAL THEORY AND APPLICATION 
TO UNRESTRICTED SYSTEMS 


By F. H. CRawrorp 
(Williams College) 
Received May 7, 1949 


I. INTRODUCTION 


The utility of Jacobians or functional determinants in obtaining 
the general first derivative relations in thermodynamics has been 
pointed out by the author! and the special case of two independent 
variables discussed in detail.2, The purpose of the present paper is to 
develop the method in terms sufficiently general to permit immediate 
application to the greatest diversity of physico-chemical systems. 
It represents essentially a unified approach to the general analytic 
treatment of thermodynamics and it is the method rather than the 
results on which any claims of novelty must rest. The method is at 
once so flexible and direct that it enables specific and precise relations 
to be written down as soon as the general expression for the differen- 
tial, dU, of the internal energy of the system can be set up and a suit- 
able set of independent variables selected. After a little practice the 
user will be enabled to write down the desired relations for any 
particular system with the minimum of intermediate equations and 
with as much confidence as he may copy them from tables in the few 
cases where they exist.® 

The general procedure is essentially an extension to the case of n- 
independent variables of the results given in Part A?. Due to the 
greatly increased complexity of the expressions involved in the case of 
n-independent variables it will be convenient to divide all the possible 
systems to which classical thermodynamics applies into two groups 
on the basis of whether there are or are not special functional re- 
strictions on certain variables which must be borne in mind throughout 


1F,. H. Crawford. Phys. Rev. 72, 521A, 1947. 

2F. H. Crawford. Am. Jour. of Physics, 17, 1-5, 1949. This will be 
referred to in the sequel as part A. 

3See P. W. Bridgman ‘‘A Condensed Collection of Thermodynamic For- 
mulas” Harvard Press, Cambridge 1925 and R. W. Geranson ‘“Thermody- 
namic Relations in multicomponent Systems” Carnegie Institution of Wash- 
ington, Washington, D. C., 1930. 
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the discussion. . All those systems where no special functional restric- 
tions exist we shall call unrestricted or simply general systems. The 
present paper will be confined to these while Part II (in preparation) 
will discuss restricted systems, 1. e., those systems characterized by 
the limited functional dependence of certain variables. 


II. NoTATION AND VARIABLES 


As in Part A where n = 2 we make our point of departure Clausius’ 
combination of the Ist and IInd Laws of Thermodynamics, 


dU = TdS + dW. (1) 


where U and S are respectively the internal energy and entropy of the 
system and dW, is the differential work performed on the system 
during the quasi-static passage along an infinitesimal length of any 
arbitrary curve of change c. For the simple system under hydrostatic 
pressure dW, = — pdV, while in general it will involve more than one 
work term, one in fact for each generalized force acting on the system. 
Each of these terms will be of the form Fd6, where F is a generalized 
force (or intensity) variable and @ the corresponding or conjugate 
geometric variable. Since TdS is of the same form we may regard 7, 
the absolute temperature, as a generalized force and S as its conjugate 
geometric variable. If now when dW, is reduced to its simplest form 


nr 
there are (n—1) work terms, i.e., dW. = & F,d,0, we have for Eqn. (1) 
j=2 


i= 


dU = TdS + Fodb, + F3d0; + ... + F,d0, (2) 
or 
n 
dU = F dé, bh Fd, -+- coe -+ F,d0, = Pe F dé; (3) 
cane 


Then Egn. (3) (or its equivalent Eqn. (2)) is the generalized Clausius 
equation and holds for any reversible differential change imaginable 
in the system. This naturally rules out of consideration non-equilib- 
rium systems as well as those in which hysteretic effects are significant. 

Since the expression for dU is assumed reduced to the smallest 
number, n, of terms compatible with the nature of the system, and 
dU is a perfect or complete differential, n is simply the number of 
independent variables involved or the total number of degrees of 
freedom of the system.* We shall designate particular independent 


4 Here n must not be confused with what we may call the “degree of varia- 
bility” as given, for example, for a chemical system by Gibbs’ phase rule. 
The latter is the number of independent intensity (or force) variables and is 
in general < n. 
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variables by x’s and any special set as (2, v2, ... 2,) and shall confine 
their choice to those in the dexter of Eqn. (3). 

A glance at this equation shows that there are 2n such variables 
(Fi, .. Fn, 0 .. 4n) and these we shall call primary variables. The 
particular primaries selected for our x-set, (x1, ... 2), is in practice 
usually dictated by experimental considerations. Thus we might 
desire to choose all 6’s (a geometric set) or all F’s (a force set) or some 
other combination. Examination soon reveals the fact that some 
types of systems permit any arlitrary choice of x’s we like while others 
donot. The first type is the unrestricted or general n-variable system. 
Such are characterized by the fact that no single primary variable 
depends on less than n others. Hence in principle at least the choice 
of an x-set is quite unrestricted. On the other hand, we have all the 
other systems where this is not true but where special restrictions on 
functional dependence must be borne in mind not only in selecting an 
a-set but in carrying out the partial differentiations concerned. 
This division is particularly desirable since systems of variable mass, 
and hence most poly-phase, poly-component chemical systems, fall in 
the restricted class and require special examination. Other simple 
examples were given in Part A (sec. VI (b)). 

There are in addition to the primary variables all the subsidiary or 
secondary quantities which are either defined explicitly in terms of the 
primaries or their differentials. The secondaries include U itself, the 
enthalpy H, Helmholtz’s free energy A, Gibbs’ function G and a 
great many others of interest in particular systems. The defining 
relations for H, A, and G, and their corresponding differentials by use 
of Eqn. (3), may be written down for reference at this point. We 
have: 


n n 
H=U-—- 2 FO and dH = TdS — 2 Od; (4) 
i= i= 
n 
A=U-—TS  anddA=— SdT + & Fid8, (5) 
i= 


n n 
G=U-TS— 2 PO and dG = — SdT — 2 dF; (6) 
i= c= 


‘It is to be noted that except for A, these are not the definitions of H and G 
which are often used, viz., U + pV and U — TS + pV, respectively. Since, 
however, all such potentials are in essence arrived at by subjecting dU in Eq. 
(3) to a Legendre transformation a very great number of such exist for the n 
variable case and there seems no particular virtue in trying to separate hydro- 
static and non-hydrostatic work terms, which the conventional definitions do. 
In any case the same set of 2n primary variables results and the applications 
of Eq. (15) in sec. V below leads to the same relations, with either choice. 
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Having selected from the 2n original primaries a set of independent 
variables (21, ... 2%n) we have n primaries left which are then the 
dependent variables. These will be indicated by (Y;... Y,). When 
we do not wish to indicate the restriction to primary variables we 
shall use (X), ... X,) for an independent set and Z; for any dependent 
variable (where 2 is now no longer confined to be n or below). Thus of 
the two general derivatives (0 Y;/0x;),, ... 2, and (0Z;/0X1)x, ... x,, 
the former may involve no secondaries, while the latter might or 
might not. Thus we have: 

The primary variables (2n in number) (Fi, Fo, .. Fn, 01, 82, . . On) 


The secondary variables U, H, A, G, ete. 
A Primary independent set (21, Yo, .. Xn) 

A Primary dependent set (*, ee 

Any independent set (X,, Xo, .. Xn) 
Any dependent set (Z,, Ze, .. Sn) 


Ill. Userut PROPERTIES OF JACOBIANS 


The ease of manipulation of Jacobian expressions depends among 
other things on certain properties among which we may list a) sign 
rules, b) reduction properties and c) transformation properties. Let 
our general n dimensional Jacobian be 


8 (Zr, ..- Zn) dZ:\ (Le aZn 
is —= (=) (=) “7 (=) = J (Zi, Zo, ee Zn) (7) 
|) a 0X, XxX, 0X, 





OZ: \\ (9Z2\... / Zn 
ar Ge ba 


OZ; . OZ, ° 
( )is short for (=) etc. and where the J notation 
X. X9 —_ Xn 


1 1 
will be used whenever the independent variables will not be in doubt. 
Here and throughout we assume that there is no functional relation 
connecting the Z’s so that J never vanishes identically. 














where 





a) Sign Rules 


Regarding J simply as an nth order determinant the following rules 
governing change of sign on altering the order of any of the Z’s 
(or X’s) follow from the well known properties of determinants. 
They are particularly useful when reducing Jacobians to some standard 
form in which the order of the Z’s and X’s is fixed. 
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Oi, .. + Be) 
The sign of J = will be 
Stia + «« Bel 

(1) changed when a neighboring pair of Z’s (or X’s) is interchanged. 

(2) changed when any pair of Z’s (or X’s) is interchanged provided 
the order of the rest is preserved. 

(3) wnchanged no matter what changes of order we make provided 
every Z appears over its original X. Thus what we may term 
“vertical pairs’ may be moved around in any arbitrary way 
without any change of sign. 





b) Reduction Properties 


Whenever a common variable occurs among the Z’s and X’s a 
reduction of the order of the Jacobian takes place. Thus in Eqn. (7) 
suppose that Z, = X;. Then the first column reduces to unity 
followed by a string of zeros. Hence on expanding by minors and 
elements of this column we have only one term, the minor after sup- 
pressing the first row and column. Hence formally we may cancel X,, 


‘1s..k&) 266... 80 
i e.. 5.5 64%... Bde 








with the proviso that X, is written outside since it must remain con- 
stant throughout. When the common elements are not a vertical 
pair they may always be cancelled with a sign alteration after inter- 
changing the appropriate Z’s (by sign rule a-2 above). 

Similarly all common pairs may be so cancelled until if there be 
(n — 1) such, we have J reduced to a single derivative. Thus 


0 CX, Xo, X3, . “a Xn; Zn) (=) (8) 

0 CX, Xo, X3 eee } oe X n) OX ,, X4X9 e * Xn—1 | 
and a partial derivative may then be regarded as simply an nth order 
Jacobian reduced by cancellation to the first order. 


This process may of course be reversed; we then “expand”’ a first 
order derivative into an nth order Jacobian (see section IV below). 








c) General Transformation Properties of Jacobians 

The chief reason for the utility of Jacobians in the present connec- 
tion lies in the ease and flexibility of their transformations. The 
Jacobian of Eqn. (7) involves n? derivatives with (Xi, X2... X,) as 
the independent set. If we wish to express this Jacobian in terms of 
a new special set (2, ... 2,) of primary variables, in terms of which 
all the Z’s and X’s are regarded as defined, we have by the well 
known theorem: 











170 CRAWFORD 


Ot, . + + Be) 





a..s0) O...4)- coh... 8) 
i...) 66...ko J... Ba 
0 (x4, -_ Za) 


Here it is quite important that we write the theorem as the ratio of 
two Jacobians—since in this case both J (Z, ... Zn) and J (X; ... Xn) 
are in terms of the same independent variable set (27 ...2,). If we 
use the mathematically identical (and usually quoted) form of the 
product of the two Jacobians, i. e., J (Z, ... Zn) XK J (a... &n), then 
the first Jacobian is in terms of the independent set (2; ... x,) while 
the second employs (X; ... X,) and the utility of the transformation 
for our application is greatly lessened. 


(9) 














IV. GENERAL DERIVATIVE AS A Ratio oF Two JACOBIANS 


We may now combine the theorem given in Eqn. (9) with the 
expansion of an arbitrary derivative, Eqn. (8), and have our general 





Z 
result. Thus, given a derivative ( ) , we express it 
2 


1 * Be 


as the ratio of two separately computable (i. e., essentially inde- 
pendent) Jacobians involving the primary set (2 ...2,) as inde- 
pendent variables. We have, expanding the derivative and using 
Eqn. (9): 

0(Z, Xo, ... Xn) 














( —) 7) (Z, Xe o 8 X») 0 (x1, wa. . - Se) 
Mije,...2, @808%...%) 00.4... 20 
) (24, ae. 2 « 2.) 


‘a. ...Ka 
’ i ie... ed 


Having selected a set of x’s we shall suppose that they are written 
in the order (21 ... 2n) in which they occur in Egn. (8) read left to 
right. Then the remaining primaries are Y’s and are to be arranged 
in a similar order (Y;... Y,). Then form the Jacobian Jo 

tts, . «+ te) 


ain Sim, .. + Mad — 


(10) 








If we call this the fundamental Jacobian for the system (its form of 
course varying with the particular 2 selection made), we have the 
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theorem: “Any arbitrary derivative ( ~-) may be ex- 
i] Xo ° xX, 

pressed as the ratio of two Jacobians which will never involve deriva- 

tives other than the n? contained in the fundamental Jacobian, Jo.” 

To prove this we shall suppose first that the (n + 1) variables in 
(0Z/0X1)x. ... x, are all primaries. Now select any n primaries as 
our x-set. Then since there are only 2n primaries all told there 
must be at least one common member in the x-set and set (Z Vo... X,,) 
or set (XY; ... X,). Hence there will be at least one cancellation in 
one of the Jacobians of (10). Suppose it is Z which cancels X, 
(it is really immaterial). Then all the derivatives in J (V2... X,)x, 
occur in J (X;, X2... X,). Since there are only n dependent varia- 
bles the set (XY; ... X,) must except for possible difference of order 
be (¥, ... Y,,). Hence no derivatives occur in Egn. (10) which are 
not in Jp. 

Now suppose any one of the set (Z, Xi, Yo ... X,) contained in our 
derivative is a secondary variable. For definiteness let it be Z again— 
then the first column of our top Jacobian will include all the deriva- 
tives of Z with respect to (7, ...2,). But since Z is expressed in 
terms of primary variables only, no new derivatives (i. e., not con- 
tained in Jo) can arise, though of course the primaries themselves may 
be involved. Hence the theorem. 


_ n(n—1) 
V. MAXWELL’s ; RELATIONS 





Thus far we have used purely mathematical transformations and 
Eqn. (3) has been used only in so far as it fixed the primary variables 
involved. We now need the most general analytical consequences of 
the thermodynamic truth contained in this equation. If we take 
this truth to be the fact that dU is a perfect or exact differential, then 
the conditions of exactness lead to precisely the type of results we 
need. But these may be written down at once by equating cross 


derivatives taking the terms in Eqn. (3) in pairs. This gives us for 
| 


OF; OF» 
example ( ) = ( ) , ete. Since there 
O05 0430, » » = @ 06; Oy + + + @ 


nm nm 








| n(n — | —P ow, 
are obviously such pairs (since the order of the pairs is im- 


n(n — 1) — 
material) we must have : such equations. This set of 


n(n — 1) 











independent equations constitutes in fact Maxwell's 
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general relations for the particular x-set (0, ... 0,). Since, however, 
we need them in the form suitable for any x-set this general form must 
now be found. This is perhaps most simply done by returning to 
Egn. (3) and expressing each dé in terms of dx; ... dz, since for each 
6,0 = 60(x,...2,). Making this substitution and collecting terms 
we have (omitting the variables held constant in each derivative) 


av - EAS Vaden + 3 Ff 


i=l i=l 








00; 
T2 

Now we apply the conditions of exactness to the coefficient of the dz’s 
taken in pairs just as before. -Thus for the first pair on differentiating 
and regrouping the uncancelled terms 


El (o.) (an) Con) Gen 


0 (F;, 6;) 


> =0 (13) 


i=1 fF) (x4, 2) 29. -* Se 














or in condensed form: 
J (Fi, 91) 25 °° * Se + J (Fo, Be) 25 oe + 
I Uatia ---2,% --- 


0 (14) 


But this states simply that Maxwell’s first relation is equivalent to the 
vanishing of n second order Jacobians in each of which (73 ... 2n) 
are held constant. Unfortunately, such a relation is too complicated 
to be readily solved for any particular derivative. We note, however, 
that the arguments in the Jacobian are simply the conjugate pairs 
(Fi, 0,) (Fe, 02), ete., as they occur in Eqn. (3). This suggests at once 
that we put a restriction on the choice of the set (a ... x,), namely, 
permit it to contain no conjugate variables, i. e., pairs of the same 
subscript such as F; and 6,, etc. The reduction ensuing is remarkable. 
In the first place J (F3, 63)2. ... 2, = 0 since x3 must either be F3 or 
6;, and every reduced Jacobian vanishes, whenever, any variable out- 
side the bracket also occurs in the numerator. Likewise J (F4, 
O4)r. ..- z, = Osince 2, is either F, or 6,, etc., for all the other terms. 
In the second place the first two Jacobians which are all that remain, 
must each reduce by cancellation to a single derivative. For since x; 
is either F; or 0; one of these cancels, and x2 being F»2 or 4 one of these 
likewise goes out. We are not able, however, to perform this cancella- 
tion until the actual explicit choice of x’s is made. Hence we must 
write as the general form of Maxwell’s first relation simply: 
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0 (Fi, 6;) 0 (Fo, 62) 


0 (x1, 2) 29. Tes 0 (x4, 2) 23. En 


= 0 (14) 








Here the Jacobians which do not vanish identically contain only 
(Fi, 0:) and (F2, 6) written in just the order they occur in Eqn. (3). 
We may then likewise write those from dx; and dz3 as 


0 (F;, 6) 0 (Fs, 63) 


6] (x4, 3) x55 rq: 0 (x4, 13) 29s ly: 





’ 





from dzz and dx; as 

















0 (Fe, 62) 0 (Fs, 63) din 
0 (22, 13) 215 4: 0 (xe, 13) 215 ry: 
or in general terms 
0 (F;, 6;) 0 (Fx, 9x) ge 1.2 ...8 
— =0 j#xk (15) 
d (x3, —~ i. 1 O (Xj, 2e)%i.-- — takes all values except j and k. 
n(n — 
The : equation contained in Eqn. (15) are Maxwell’s rela- 
tions in their most general form for the case of (a2; ...2,) being 


non-conjugate primary variables.’ In all cases the relations of Maxwell 

connect elements on one side of the positive diagonal of Jo with their 

conjugate elements on the other side, in fact the conjugates are always 

either equal in pairs or differ only in sign.” Consequently of the n? 
- _ n(n — 1) 

original derivatives or elements in Jo, : may be eliminated, 


n(n — 1) n(n + 1) 
leaving us n* — : or 
terms of which any other derivative can be expressed by means of 


Eqn. (10). Of course in many cases there will be a good deal of 


n(n + 1) 





as the irreducible number in 











reduction in the two Jacobians involved—but no more than 
need ever occur. Hence the final result: “Any given derivative, 


6 This result is quite similar to the one given for n = 2, Part A, sec. III (b) 
but the author has been unable to deduce it by such a concise method. 

7 The term conjugate is used in two senses, one referring to conjugate 
variables in Eqn. (3), i. e., variables from the same term Fjd6; and conjugate 
elements in a determinate, i. e., element pairs which are symmetrically placed 
with respect to the main diagonal. 
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OZ 
(=) can be expressed as the rat of two separately com- 
|S fa oe 

n(n + 1) 


putable Jacohians which will not involve more than : members 





of a standard set of derivatives.” 
It is clear that this standard set of derivatives must always contain 
n(n — 1) 





the n elements on the positive diagonal of J) together with 


off-diagonal elements (and not, of course, including any conjugate 
pairs). Just which are selected for the standerd set depends on con- 
venience (usually experimental!) though the terms in Eqn. (3) may 
often be arranged so as to have the standard set on and above the 


n(n + 1) 


ee 





positive diagonal. It will be noted that there are precisely 


elements in such a triangular array and once the elements of such an 
array are determined experimentally, all other derivatives are in 
theory fixed. Since this triangular array is really the functional 
Pfaffian of the dependent with respect to the independent variables, 
we may for brevity refer to this array simply as the Pfaffian of Jo. 
No reference is intended to the actual magnitude of the Pfaffian as an 
algebraic quantity, since in our reckoning square arrays only may 
arise. On this basis restricted systems will be seen to differ from the 
general unrestricted ones chiefly in the fact that for them one or 
more elements of the Pfaffian of Jo must vanish. 
VI. Summary oF METHOD FOR UNRESTRICTED SYSTEMS 
Although the above treatment covers the unrestricted system 
where primary variables offer suitahle x-sets, there are a number of 
points which in practice greatly simplify the whole discussion. 
The simplest case arises when all the terms in Egn. (3) are positive 
and we choose the x set as all F’s or all @’s._ In these cases Eqns. (15) 


become, respectively, 
06; 06}. 
( , ) _ ( ) (16) 
OF k F. OF, F. 








and 








OF; = OF, i1# k . 
00, Jo. 00; J 6. 2 includes all integers but 7 and k (17) 


But this means that, with either choice, Jo will be symmetric about the 
main diagonal. To be specific let the x’s be all forces (when writing 
S for 0, and 7 for F;) we have for J, 











— 
~J 
qn 
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P 0 (S, Oo, . . . On) 
v3 = 
O(T, Fo, ... F,) 
OS 06. 06s 06, 
o dg@ |— — —...—| = | @ @ ..-. Gmi was 
OT oT oT oT 
OS 00. 063 
- oe GQ, @j,; Gyo ... Aim 
OF, OF, OF» 
OS 062 063 
ai qlee i (lo Qjo Gog . . « Qom 
OF; OF; OF; 
OS 06, 
ny (2a | ics. «a os Sie 
OF, OF,,| 
Here the two square arrays one of derivatives and the other of coeffi- 
cients do... Anm (with m = n — 1) serve as a convenient tabular dis- 


play giving us Maxwell’s relations at a glance as well as a concise set 
of abbreviations for the derivatives proper.* The derivatives nat- 
urally break up into two groups, the first column being calorimetric 
derivatives and the rest of the first row thermal dilations. All of the 
former except (0S/0T),° are equal by symmetry to corresponding 
members of the latter group; since the thermal dilations are usually 
more readily measured they are usually retained and the coefficients 
in Ao, @1 ... Gm, would accordingly be taken as representing gener- 
alized thermal dilations. The remaining square array of m? isothermal 
derivatives are reduced by symmetry to m (m-+ 1)/2 independent 
generalized isothermal elastic strain coefficients. Just which of an 
equal pair is to be eliminated in practice depends on circumstances 
and the method of writing Ap leaves the choice of course open. 

As we pass from systems with small values of n to more complicated 
ones a new and longer column is added to the Pfaffian for each unit 
increase in n, as the table below shows. In the first column is the 
value of n, in the second the number of independent coefficients in the 
fundamental set (or in the Pfaffian of Jo) and in the third a typical 
example. Of these the case of n = 2 has been considered in Part A 
and for the last three it will be sufficient for our purposes to set up Jo 
or its equivalent as an indication of both the closed nature of the 
analysis and of the way in which the results of more specialized 
treatments are summarized in the general case. 


§ When it is necessary to distinguish between the two square arrays as such, 
Jo and Apo will be used. 

*Since this is (0S/08T)F, . . . Fn, @ is the heat capacity of the system at 
constant F, ... Fn divided by the absolute temperature. 
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TABLE 
Derivatives 
in 

n Pfaffian Typical System (unrestricted) 

2 3 Pure substance under hydrostatic pressure 

3 6 Uniform Rod under Tension and pressure 

4 10 Triclinic crystal in a uniform Electric field 

7 28 Triclinic crystal under uniform general stress 
10 55 Triclinic crystal under stress and Electric field 





It is of course an advantage to have the fundamental Jacobian 
symmetric since this means that we may write down Maxwell’s rela- 
tions directly by equating conjugate derivatives without using Egns. 
(15,) (16) or (17) at all. There are two remaining cases, which in- 
clude the most important practical situations, where the symmetry 
may be retained by a simple rule of signs. ‘Thus whenever any term in 
Eqn. (3) appears with a minus sign we insert a minus sign in front of 
the independent variable in Jp. (The most common case of course is 
that of the hydrostatic term occurring as — pdV.) Further, if for 
some reason we must depart from all-force or all-geometric variables, 
i. e., use a “mixed” set of primaries, we may preserve our symmetry 
by introducing a minus sign in front of any independent variable 
(in Jo) whose order in the terms Fdé@ is opposite that of the order 
established in selecting from F; and 6; in F;d6,._ Thus if we take for 
our x-set the first variable in the first term of Eqn. (3), i. e., Fi, but 
select the second (62) in the second, F2 d@2, we must put a minus sign 
in front of 6 in Jo.!° Obviously, if we are using all force variables 
except in a hydrostatic term where we choose V instead of p, the two 
sign changes cancel and + V is used in Jp. These sign rules will 
become somewhat clearer after the first example below has been 
discussed. 

Before beginning the examples, one further remark should be made. 
The whole of the discussion so far could have, with equal ease, been 
based on the equation for dH, dA or dG, Eqns. (4), (5) or (6). Thus 
each gives the same set of primaries for the problem. Further, if we 
apply Maxwell’s relations in the general form of Eqn. (15) we obtain 
precisely the same results, provided that the terms F; d0; or 0; dF; are 
written in the same order as in Egn. (3). This results from the fact 
that J (F;, 6;) = — J (0;, F;) and that the terms in the four Eqns. 


10 Tt really is not important whether the minus sign occurs before the de- 
pendent or the independent variable since in either case we change the sign of 
all the elements in a row or column and the formal symmetry of the Jacobian 
is preserved. 
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(3) to (6) are in fact always identical except when force and geo- 
metric variables are interchanged (when minus signs always occur). 


VII. EXAMPLeEs 


Of the unrestricted systems met in practice, probably the most 
common is that for n = 3, and it will be treated at slightly more 
length than the others. 


a) n = 3 Thin rod under pressure and tension 


Given a thin rod of volume V and length L which is under both 
uniform hydrostatic pressure p and an extensive force F (gravitational 
effects being neglected). Eqn. (3) becomes 


dU = TdS — pdV + FaL (19) 


where we have the 6 primary variables (7, S, p, V, F, L) any three of 
which we may in theory take as independent. If we choose (7, p, F) 
for our x-set we should then write 


0 (S, V, L) os OV ol 

Jo = = omnes aemmnene —— | = | ao a ae = Ao 
a(T, — vp, F) aT oT aT 

os OV OV | 

Sa a ee ee aa A, Ay, Aye (20) 

Op Op ap 

os OV OV 


OF oF oF 








Qa Ajo Ax 














where the (— p) is used to preserve the symmetry of J) and give us 
the reciprocal relations directly on equating conjugate derivatives 


OV os ian 
(as (=) =— (—-) , ete.) Here the a’s are abbreviations for the 
oT pF Op/ rr 


corresponding derivatives and enable us to eliminate the three 
unnecessary derivatives (usually those below the main diagonal in J») 
in any particular calculation. On the other hand, in general calcula- 
tions it is awkward to have to keep inserting this minus sign in front 
of p whenever it occurs. We may accordingly avoid this trouble by 
taking a new Jacobian, Jp’ = 0 (S, V, L)/d (T, p, F), as our funda- 
mental one. On writing Jo’ out we set it equal to a symmetric array 
of a’s (Ap in Eqn. (20)) and then insert minus signs according to the 
two sign rules in Ap only (here in front of the a’s in the middle row. 
which correspond to the p derivatives). Thus we have 
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0d (S, V, L) OS OV ab 
Je = — = — —/={| © Q om] =A 

0 (T, p, F) OT OT oT 
OS OV aL 
ee. ee dé, ~Ga, ~~ Gan (21) 
Op Op Op 
OS OV aL 
— | & a a 
‘eaii"” ™ 











In all cases Jo’ = + Jo (according to the number of sign changes 
introduced), though the important fact is that now Ao’ is always a 
symmetric determinant while Jo’ may or may not be. 

It can now be shown that every derivative made up of primaries 
only, can be expressed as a ratio of Jacobians which can be read 
directly out of Jo’ (so that tables of evaluated Jacobians need never 
be prepared unless secondaries are in frequent use). Thus the only 
Jacobians which can arise are those of the form (J (11, Yo, Y3), where 
the Y’s are to be selected in all possible ways from the primaries, 
(T, S, p, V, F, L). But there are only 6.5.4/3! = 20 essentially 
different ways of doing this (since we divide by 3! to eliminate the 
result of permuting the order of the Y’s). There are thus 20 Jaco- 
bians all told; one of those is J (7, P, F), (which is unity), 9 are the 
elements in Jo’, 9 are the various minors formed by suppressing a 
single row and column and the last is Jo’ itself. 

As an example suppose we wish the heat capacity of the system at 
constant L and V. Since Cry = T (dS/dT) zy we proceed to obtain 
this derivative. On expanding the derivative into a third order 
Jacobian and transforming as in Eqn. (10), from (7, V, Z) as inde- 
pendent variables to the standard set (7, p, F) we have: 


0 (S, V, L) 


as as 0d (S, V, L) d(T, p, F) 
aa). 7 Gal. a7, V,L) a(T,V,L) 
d(T, p, F) 
Jo +J0 


J’ (T, V, L) J'(V, L)r 
But Jo’ we have and J’ (V, L)r is read from Jo’ by suppressing the 
S-column (since S is not involved) and the 7-row (since T is constant). 
Therefore 




















Ay A, Ae | 
Q Qj, Qy2! 
Ag Ajo Age! 








Ciy = Crit = 


en Qy2! 
Ayo Ago! 


where the two minus signs have cancelled. 
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Again to find the effect on temperature of adiabatic stretching at 
constant p we want (07/0F),s. Proceeding as before 

















d(T, p, S) 

oT d(T, p, S) d(T, p, S) d(T, p, F) 
Ga). “sana ~ 8(p,F,8) —-a(p, F,8) 

d(T, p, F) 


as 
(oF), 

as 
(or), 





and again reading directly from Jo’ in Eqn. (21) 


(=) ay"! 
oF ps Ap ) 


The ease and economy of such deductions are typical of the method 
and are in marked contrast with deductions by conventional methods, 
as may be seen by a trial. 

If for some reason or other we wish to select another z-set such as 
say (7, p, L) we have 











a(S, V, F) aS OV OaF 
J; = = —— a sae == bo by be = A)’ 
d(T, p, L) oT oT aT 
aS aV OF 
oe eee —b, —by — bie 
Op Op op 
aS aV oF a 
snl psd es 
aL aL aL ee 














where 6’s are used instead of a’s since these derivatives are different 
from those in Eqns. (20) and (21) and the minus sign is inserted in 
the middle row due to the sign of pdV and in the bottom row since L 
was taken out of turn from the term FdL. Use of Eqn. (15) verifies 


av as as 
the fact that indeed (=) = — () ( ) on 
oT pL Op /TL OL Tp 


OF 
—{_— , etc. 
Gs) 


11 This equation then tells us that if F be increased the change in T is op- 
posite in sign to az (since a» is surely positive) and therefore normal substances 
cool while materials such as stressed rubber, where az < 0, warm on stretching. 
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Again if we chose (7, V, F) we should have: 





d (S, p, L) OS dp OL 
J, = = —— aE ones te tia oe = A,’ 
d(T, V, F) OT OT OT 
0S dp OL 
aV aV av a 
OS dp dob 
— —— —— c 
aF OF dF aliens 














where here both Jo’ and Ay’ are symmetric. The use of such mixed 
z-sets usually arises in restricted systems where for example F or p 
may be a function of 7’ only and hence inappropriate in an z-set con- 
taining 7 itself. In this connection it should be emphasized further 
that systems will arise in practice where only one (usually 7) of the 
primary variables is convenient for inclusion in the x-set. In this 
case the above arguments based on Eqns. (15) ceases to apply and we 
must go back to the general Maxwell’s Eqns. as given for example in 
Eqn. (14). It can similarly be seen that Eqn. (14) is still valid even 
when the expression for dU has not been reduced to the minimum 
number, n, of differential terms. If for example there are n’ (n’ > n) 
terms in Egn. (3) we have n’ second order Jacobians in Eqn. (14) and 
although the result looks cumbersome it is nevertheless more con- 
venient in practice than the usual procedure in which one essentially 
has to repeat the whole process of Section V in each individual case. 
Where secondaries are involved the Jacobians concerned will be 
more complicated in that derivatives of secondaries with respect to 
(T, p, F) (or whatever x-set is being used) will appear. These may 
readily be got from Eqn. (19) and the corresponding forms of Eqns. 
(4), (5) and (6). “‘Q” and “W’” mav be allowed among the secondaries 
formally since dQ. = TdS and dW, = — pdV + FdlL. We have 
simply to replace any partial derivative which occurs formally in a 


dq 
Jacobian by its equivalent. Thus for example (=) is interpret- 
P /TV 
ed as dQry/dp, where this symbol represents the derivative of the 
heat function at constant 7 and V and may then be replaced by 


T (0S/dp) rv. 
b) n = 4 Single triclinic crystal in an electric field (p = 0). 


Consider a single crystal (of lowest symmetry) in an electric field 
of components EF), E2 and E; (these are values of the electric field in- 
side the crystal). Let the corresponding components of the electric 
induction be D,, D. and D3. Then the expression for dU, the differ- 
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ential of the internal energy of the system per unit volume, can be 
written in the form!” 


Strictly speaking we should then divide by the mass density, p, and 
consider du = 1/p dU, where u refers to a fixed mass of unity, as the 
analogue of Eqn. (3). Actually, since the volume changes due to 
electrostriction are quite small, Eqn. (22) is adequate. With n = 4 
the Pfaffian contains 4.5/2 = 10 terms of which 4 are thermal and 6 
isothermal, as the following considerations show. 

Since in such a general case the D vector will not be parallel to the 
E vector, each component of D will have contributions from F,, FE, 
and E3. Hence 


D, = Ky Fi + Ky ke + Ky; Es 
Dy = Ko Ey + Koo Fo + Kop Eg. . (23) 
D3 = K3, EF; + Ky Ex + Kzs Ez 


where the K’s are the dielectric constants of the crystal. 
On writing Jp we have 


























d (S, Dy, De, Ds) 0S OD, IAD. ADs 
0 = = a ae (23A) 
0 (T, Ey, Es, E3) OT OT oT oT 
oS oD, 
OF, OF, 
oS oD, 
OE, OE, 
oS aD, 
OF; OF; 








But the symmetry requirements reduce to K;, = Ky; for the K’s in 
Eqn. (23) and we have 


Jo= |Q9a% Mm a3 | = Ay 

ay Ku Ky Ki e ‘ ; (24) 
dz Kio Kee Kos 
dg Ky3 Kos Ka 








12See, for example, JEANs: “Electricity and Magnetism.” Cam. Univ. 
Press, 4th ed., page 152, where a slightly different notation is used. Unless 
the crystal in this and the examples below belongs to the lowest class of 
symmetry, the various symmetry demands of the other classes change the 
system from an unrestricted to a restricted one and therefore modify the 
treatment considerably. 
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| — . , 9D: aD: 
Here as usual ay is a calorimetric coefficient, while ar » —— ,ete., 
0 


oT 


are pyro-electric coefficients, and the isothermal coefficients are simply 
the 6 independent dielectric constants of the triclinic crystal. 

Now it is seen at once that the three pyroelectric coefficients vanish 
with the field. Thus from Egn. (23) to (24) we have 


OT Je,..-n, \OT aT aT 


(=r) 
OF, /re.e, 


and hence there is no pyroelectric effect unless there is a field present. 
Hence no pyroelectric effect in the usual sense of the word is present 
unless we have added to each of the right sides of Eqns. (23) terms 
Dy, Dey and Dz, respectively, where these quantities represent 
permanent electric moments (being functions of temperature in this 
case—but of stress as well in the general case). These permanent 
states of polarization are then the origins of true pyroelectricity and 
can be expected only for crystals having a onewayness in some direc- 
tion, i. e., at least one polar axis due to the physical asymmetry of the 
crystal class. This is of course well known, but the above analysis 
gives the fundamental thermodynamic requirements in a very direct 
way. 








c) n = 7 Triclinic crystal under small homogeneous strain 


If we now imagine our single crystal to be subject to a general 
stress it can be shown" that Eqn. (3) for a unit volume becomes 


6 
dU = TdS + >° F;d8, (25) 
i=1 


where (F, ... Fs) are the usual components of stress (X:, Xy, V2, Z:, 
X,) while (6; ... 65) are the corresponding strain components (ézz, 
Cyyy Cee) Cys) Coz ANd Czy). 


In this case the generalized Hooke’s law enables us to express the 
6’s directly in terms of the F’s. We have in fact 6 equations (analo- 
gous to the 4 in Eqns. (23)) expressing each @ as a linear function of 
the F’s. Thus 


13 See A. E. H. Love, Mathematical Theory of Elasticity. Camb. U. Press, 
4th ed., 1934, page 94. 
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6, = 00 + Cy Fy + Cp eo + Cig '3 + Cu Fs + Cis Fs + Cre Fe 
(26) 


62 = Oo + Ca Fi + Cox Fo + Coz Fs + Cog Fa + Cos Fs + Cog Fo, ete. 


where 640, 40, etc., are functions of temperature only. These involve 
36 general elastic constants Cy... Ce. Just as before we set up Jo 
which is a square array of 7? = 49 terms consisting of a smaller square 
of 36 2sothermal constants bordered by a row and column of 13 thermal 
and entropy constants. By the symmetry requirements these reduce 
to 7.8/2 = 28 independent terms, 7 thermal and 6.7/2 = 21 isothermal. 
Thus putting Cj, = C;;, etc., we have 

0(S, 0. . . 4) Ma & 


o © « & 
= = |aCuCn. . . 
O(T,Fi--- Fe) laCnCn. 2 | Cm) 2D 





0 








ae 
As before we are now in a position to determine in terms of these 28 
coefficients any new derivative. As an example consider any dy- 
namically determined elastic constant where the deformations are so 
rapid that we may assume them adzabatic, such as (062/0F¢) SFF3- + - Fe 


Then 
00, _ 8 (6, S, Fi, Fs. . . Fe) 


Ge) nem MAA... .0 
0 (S, 02, Fi, Fs. . . Fe) 
0 (S, Fi, Fo, Fs . . . Fe) 
0 (S, 6, Fi, Fs. . . Fe) 
Ay Ap | 
0(T, Fo, Fi, Fs. . . Fe) J(S, 02) ry mer 3 &% Ae Cr 
~ 8(S8, Fi, Fo . ‘ . Fe) (=) A 
0(7, Fi, Fo... Fe) F, 























oT - Fe 
where a number of short cuts suggest themselves automatically in 
practice. 


d) n = 10 General Piezo electric crystal 


It is apparent without further repetition of similar details that we 
have the material at hand for a complete formal discussion of the 
general piezo electric crystal, i. e., of the triclinic crystal under general 
mechanical stress and arbitrary electric field. If we add to TdS the 
non-thermal parts of Eqns. (22) and (25) we have for a unit volume 
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6 3 
qU = TaS +E F,do, +E EdD; (28) 
{= i= 


Since n is now 10 the Pfaffian contains 10.11/2 = 55 coefficients and 
without writing them‘out in detail their general nature becomes clear 
from the block diagram of the 10th order Pfaffian. On the top row 





—" 7 3 





21 18 








we have a and 6 dilations from Eqn. (27) and the 3 pyro-electric 
coefficients from Eqn. (23A). Below are triangular blocks of the 21 
isothermal elastic constants and 6 dielectric constants from the same 
equations and finally a rectangular block of 18 new coefficients the 
piezo-electric constants proper. ‘These latter result from the inter- 
action of the stress and the field and bring the total to the number 
required by our theory, viz., 55." 


14 See for example W. G. Cady: ‘‘Piezoelectricity.”” McGraw-Hill, 1946, p. 
44 where essentially the same group of constants is given (in quite a different 
notation) in their relation to a closed expression for A, the work function. 











